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Distinguishing the nontrivial symmetry-protected topological (SPT) phase from the trivial insu¬ 
lator phase in the presence of electron-electron interaction is an urgent question to the study of 
topological insulators, due to the fact that most of the topological indices defined for free electron 
systems are very likely unsuitable for interacting cases. In this work, we demonstrate that the 
strange correlator is a sensitive diagnosis to detect SPT states in interacting systems. Employing 
large-scale quantum Monte Carlo (QMC) simulations, we investigate the interaction-driven quan¬ 
tum phase transition in the Kane-Mele-Hubbard model. The transition from the quantum spin 
Hall insulator at weak interaction to an antiferromagnetic Mott insulator at strong interaction can 
be readily detected by the momentum space behavior of the strange correlator in single-particle, 
spin, and pairing sectors. The interaction effects on the symmetry-protected edge states in various 
sectors, i.e., the helical Luttinger liquid behavior, are well captured in the QMC measurements of 
strange correlators. Moreover, we demonstrate that the strange correlator is technically easier to 
implement in QMC and more robust in performance than other proposed numerical diagnoses for 
interacting topological states, as only static correlations are needed. The attempt in this work paves 
the way for using the strange correlator to study interaction-driven topological phase transitions in 
fermionic as well as bosonic systems. 

PACS numbers: 71.10.Fd, 71.27.+a, 73.43.-f 


I. INTRODUCTION 


Topological insulators (TIs) are usually defined as sys¬ 
tems with bulk spectra similar to those of trivial in¬ 
sulators but with nontrivial, i.e. gapless or degenerate, 
boundary spectrums when and only when the systems 
(including the boundaries) preserve certain symmetries. 
By now the noninteracting TIs have been fully classi¬ 
fied and understood, for example, as shown in Refs. 1 
3. Besides the boundary states, which are experimen¬ 
tally most relevant, the noninteracting TIs can also be 
characterized by a topological index defined for the bulk 
band structure, namely, even if a TI has the similar bulk 
spectrum as a trivial insulator, it does have a very dif¬ 
ferent ground-state wave function which is characterized 
by the topological indices, for example, the Thouless- 
CKohmotoCNightingaleCden Nijs (TKNN) number for 
the integer quantum Hall state 4 and the Z 2 index for 
the quantum spin Hall insulator 5,6 . To generalize the 
notion of topological insulator to interacting systems, 
symmetry-protected topological (SPT) order 7,8 was pro¬ 
posed for the ground states of many-body quantum sys¬ 
tems that have a symmetry and a finite energy gap with 
short range quantum entanglement. So far most of the 
techniques and topological indices introduced for non¬ 
interacting topological insulators are very likely unsuit¬ 
able for interacting cases, since in many cases interaction 
can change (or reduce) the classification of topological 
insulators 9-19 . Thus a more general technique to iden¬ 
tify interacting TIs (or SPT states) based on their bulk 
wave functions is urgently demanded for studying topo¬ 


logical insulators in the presence of electron-electron in¬ 
teractions. 

In principle, given a bulk wave function, we can always 
compute its entanglement spectrum and use it as a di¬ 
agnosis for interacting TI 211 . However, this technique is 
numerically challenging. For strongly correlated electron 
systems, in one dimension (ID) we are able to obtain the 
bulk wave function and entanglement spectrum from ex¬ 
act diagonalization (ED) and density matrix renormaliza¬ 
tion group (DMRG) calculations, but in two dimensions 
(2D) and higher, it is very difficult to obtain the bulk 
wave functions for interacting systems simply because 
the dimension of Hilbert space increases exponentially 
with the number of electrons. In 2D, there has been re¬ 
cent progress by employing quantum Monte Carlo simu¬ 
lations to access the entanglement spectrum 21-23 , but the 
approach is arduous, as one needs to first bifurcate the 
already small finite-size system (the simulation efforts of 
QMC scale polynomially with system size to high power) 
and then perform analytical continuation to obtain the 
real-frequency entanglement spectrum from the reduced 
density matrix in imaginary time 22,23 . The analytical 
continuation 24,25 , as useful as it is, is a numerically ill- 
posed question and is used with caution for bringing am¬ 
biguities that mask the fine features in the real frequency 
data. These difficulties shadow the progress in evaluating 
the bulk wave functions and entanglement information 
for diagnosing interacting TIs. 

In light of the difficult situation for interacting TIs, 
recently, a new diagnosis dubbed “strange correlator” 
was proposed in Ref. 26, which is the matrix element 
of the correlation function between two topologically dis- 
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tinct many-body bulk wave functions in the same Hilbert 
space. Based on the low-energy effective Lorentz invari¬ 
ance of the SPT states, the strange correlator, though it 
is a purely static quantity, effectively captures the space- 
time correlation function at the spatial interface between 
two topologically distinct phases. Hence, as long as there 
exist symmetry-protected edge states at the spatial in¬ 
terface, i.e., the two wave functions are topologically dis¬ 
tinct, the strange correlator will diagnose the edge modes, 
at least for the noninteracting case. As will become clear 
in this paper, for the QMC simulations of interacting 
TIs, the strange correlator can diagnose the correlated 
edge modes as well. Moreover, as the strange correla¬ 
tor is based on the bulk wave function, there is no need 
to explicitly create a real spatial boundary to detect the 
gapless edge modes, which, in interacting systems, usu¬ 
ally gives rise to strong finite-size effects. Also, com¬ 
paring with the measurements of the entanglement spec¬ 
trum mentioned above, there is no need to bifurcate the 
system for evaluating the strange correlators. There is 
also no need to perform imaginary-time correlation as the 
strange correlators are static quantities which avoids the 
analytical continuation step. These advantages make the 
strange correlator physically transparent and technically 
much easier to implement in QMC. Yet another advan¬ 
tage of the strange correlator is that it is generally appli¬ 
cable to both fermionic and bosonic SPT states, either 
free or interacting. It is also applicable to “crystalline” 
TI 27 , because it respects all the lattice symmetries (no 
need for boundaries). 

In Ref. 26, the strange correlator has been applied only 
to free fermion topological insulators and some bosonic 
SPT states. Later on, in Ref. 28 and 29, it was demon¬ 
strated that the strange correlator can capture the nature 
of the Haldane phase of ID spin-1 systems. It was fur¬ 
ther shown in Ref. 30 that the strange correlators of 2D 
bosonic SPT states can be expressed as correlation func¬ 
tions of 2D conformal field theory. However, the most im¬ 
portant test, namely, the application of the strange cor¬ 
relator upon interacting fermion topological insulators to 
diagnose the interaction-driven topological phase transi¬ 
tion, has never been performed. Here, by means of large- 
scale quantum Monte Carlo simulations, we apply the 
strange correlator to a very realistic and nonintegrable 
model for interacting topological insulators, namely, the 
Kane-Mele-Hubbard (KMH) model. We present details 
on how to evaluate the strange correlators in determi- 
nantal QMC 31 simulations for interacting fermionic sys¬ 
tems and use it to probe the topological nature of the 
interaction-driven quantum phase transition in the KMH 
model. Furthermore, the interaction effects on the heli¬ 
cal edge states - the Luttinger liquid behavior have been 
also clearly captured by the strange correlator measure¬ 
ments in QMC simulations. 

The rest of the paper is organized as follows. In Sec. II 
the KMH model (II A) and strange correlators in various 
sectors (IIB) are introduced, with detailed accounts of 
their implementation in QMC simulations. In Sec. Ill, 


the strange correlator in the single-particle sector (III A) 
is first demonstrated, followed by those in two-particle 
spin and pairing sectors (IIIB). In the single-particle sec¬ 
tor, the topological nature of the quantum spin Hall insu¬ 
lator to antiferromagnetic Mott insulator transition can 
be clearly seen. In the two-particle channel, the Lut¬ 
tinger liquid behavior of the edge modes, is well captured 
by their corresponding QMC strange correlator measure¬ 
ments. Section. IV summarizes the physical and numer¬ 
ical advantages of the strange correlator in diagnosing 
interacting TIs and proposes future directions. 

II. MODEL AND NUMERICAL METHOD 
A. generalized Kane-Mele-Hubbard model 

The generalized KMH model is given by, 

^KMH = tij c icr C jcr + *A v ij c \ a a a.p c iP 

+IB"*- 1 ) 2 - a) 

i 

Here the first term describes the nearest-neighbor hop¬ 
ping on a honeycomb lattice. As shown in Fig. 1 (a), we 
set the nearest-neighbor hopping within one unit cell with 
amplitude td, while others are set with amplitude t, t^ 
and t can be different, depending on the context. The sec¬ 
ond term represents spin-orbit coupling (SOC) 5 ’ 6 , which 
connects the next-nearest-neighbor sites with a complex 
(time-reversal symmetric) hopping with amplitude A, and 
the factor = — Vji = ±1 depends on the orientation 
of the two nearest-neighbor bonds that the electron tra¬ 
verses in going from site j to i. The cC 3 in the spin-orbit 
coupling term furthermore distinguishes the f and j, spin 
states with the opposite next-nearest-neighbor hopping 
amplitudes. 

Physically, the noninteracting (U = 0) Kane-Mele 
(KM) model 5,6 is a spinful model consisting of two copies 
of the the Haldane model 32 with opposite spins. Al¬ 
though the spinless Haldane model breaks the time- 
reversal symmetry Zj, the spinful KM model is time- 
reversal invariant and its ground state is a quantum 
spin Hall (QSH) insulator with counter propagating edge 
modes. Regarding the symmetries of the KMH model, 
the model Hamiltonian in Eq. 1 has the charge U{ l) c harge 
symmetry Cj CT —> e lS Ci a . The spin-rotational symmetry 
SU( 2) is broken down by the spin-orbit coupling term 
A to U (l)spi n , which keeps the spin rotation only in the 
xy plane: Ci a e lcre Ci a . So the symmetry group of the 
KMH model is U(l) c harge x C7 (l) sp in x Zj, which results 
in a Z classification. This implies that the QSH state 
must be separated from the trivial vacuum state with 
gapless edge modes. 

In the presence of interaction, the KMH model can 
be studied by determinantal QMC simulations 31,33-40 . 
At the noninteracting limit, U = 0, for any finite A, 
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FIG. 1. (Color online) (a) Illustration of honeycomb lattice 
and Kane-Mele model. The unit cell of the honeycomb lat¬ 
tice is presented as the dashed black parallelogram, it consists 
of two sublattices, A and B, denoted by the open and filled 
cyan circles. The underlying lattice is spanned by the prim¬ 
itive vectors ai = (%/3, 0), <22 = (-s/3/2, 3/2). The green and 
orange lines represent nearest-neighbor hopping t and td con¬ 
necting the A and B sublattices. The spin-orbital coupling 
term (complex valued next-nearest-neighbor spin-dependent 
hopping iX) connects lattices sites within the same sublattice 
and is denoted as the blue dashed arrows, (b) Schematic plot 
of |Q) and I'F) in single-particle strange correlator calculation, 
where |fi) is a trivial band insulator with spin Chern number 
C s = 0, and the detected target state |'F) is the many-body 
ground wave function of the KMH Hamiltonian in Eq. 1 eval¬ 
uated in the QMC simulation. It is prepared by applying the 
projection operator e~ eH onto a noninteracting trivial wave 
function |<E't) (eigenstate of the KM model with U = 0). The 
topological nature of |$) depends on the interaction strength 
U/t. When U < U c the system is in the QSH insulator phase 
with spin Chern number C a = 1, but when U > U c , the sys¬ 
tem is the AFMI phase, which spontaneously breaks the key 
symmetry that protects the topological insulator. 


the system is in the QSH state at zero temperature. 
Switching on finite but weak interaction (U/t > 0), 
the system is adiabatically connected to the noninter¬ 
acting QSH. At strong interaction, U/t will drive the 
QSH state into an antiferromagnetic ordered Mott in¬ 
sulator (AFMI) state 41 through a continuous quantum 
phase transition at critical point U c (e.g., at A = O.lf, 
U c ~ 5 1). At the transition, the single-particle gap re¬ 
mains open but the corresponding spin gap closes 35,37,40 . 
The transition from the QSH to the xy AFMI has been 
shown to be consistent with the 3D XY universality 
class 35,36,42 . As both t7(l) sp i n and Z J symmetries are 



FIG. 2. (Color online) Illustration of the k mesh in the Bril- 
louin zone (BZ) of finite-size systems studied in QMC sim¬ 
ulations with linear system size L = 9,12,15,18. The black 
line is the high-symmetry paths F — > M —> K — >■ F. As 
L increases, the k mesh becomes denser. The inset is the 
hexagon BZ of the honeycomb lattice, and the shaded region 
represents the segment of the BZ shown in the main panel. 


spontaneously broken in the AFMI phase, only f/(l) c har g e 
remains. Meanwhile there is another time-reversal-like 
symmetry Zj : Cj a —> JCa^pCip, so the total remaining 

symmetry is 17(1)charge x Zj with T 2 = 1. Thus the 
fermion SPT classification becomes trivial, and as such, 
a time-reversal-like symmetry with Y 2 = 1 does not lead 
to a Kramers doublet. This means that if we neglect the 
Goldstone mode in the bosonic sector, the AFMI state 
must belong to the trivial SPT class in the fermionic sec¬ 
tor, which can be smoothly connected to a trivial band 
insulator [such as a spin density wave (SDW) insulator]. 
So there is no symmetry-protected gapless fermionic edge 
mode between this AFMI and a trivial insulator. 

Of course, interaction-driven topological phase tran¬ 
sitions happen in other models as well. For example, 
in the Bernevig-Hughes-Zhang model, dynamical mean- 
field theory studies 43-46 reveal interesting (first-order) 
topological phase transitions in the paramagnetic sector 
of the solution. 


B. Strange correlator in QMC 


To effectively diagnose the SPT states, the concept 
of a strange correlator was proposed in Ref. 26. It is 
a correlation function defined between two many-body 
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wave functions in the same Hilbert space, 


C(r,r') 


{®\<t>{r)<l>(r') |T) 

<w 


( 2 ) 


where |fi) is a trivial band insulator, and |T) is the wave 
function whose topological nature we would like to di¬ 
agnose. The physical meaning of C(r, r') becomes man¬ 
ifest after a space-time rotation 26 : |T) can be obtained 
by evolving a generic initial state from imaginary time 
r = -oo to r = 0 with the parent Hamiltonian of |T), 
and |fi) can be obtained by evolving a generic final state 
backward from imaginary time r = +oo to r = 0, thus 
C (r, r’) can be viewed as a correlation function at the 
temporal domain wall. Because most of the topologi¬ 
cal insulators have an effective Lorentz invariant descrip¬ 
tion 47 , after a space-time rotation C(r,r') becomes the 
space-time correlation at the spatial interface between 
14’} and |fi), which may have gapless modes depending 
on the nature of these two states. 

The proposition given in Ref. 26 is that if IT) is a non¬ 
trivial topological insulator (or more generally, a SPT 
state) in one or two spatial dimensions, i.e., there exists 
one or more gapless edge modes at the spatial bound¬ 
ary of |T), then for local operator cf>(r) that transforms 
nontrivially under symmetry, C(r,r') will either develop 
long-range order (saturate to a constant ) or decay as a 
power law in the limit |r — r'\ —> + 00 , which mimics the 
edge states of |T). In the momentum space, this corre¬ 
sponds to a singularity at a certain symmetric momen¬ 
tum point k s : Ck ~ 1/1 k — k s |“, if |T) is in a nontrivial 
topological insulator phase. Based on the space-time ro¬ 
tation argument given above, the 2D strange correlator 
C(r,r') should behave very similarly to the (1+1)D cor¬ 
relation functions at the boundary. For example, if |T) 
is a generic noninteracting 2D TI, and cf)(r) is simply the 
electron operator, i.e. C(r,r') = (filet (r)c(r , )|T)/(fi|'F), 
then a = 1. The strange correlator has been successfully 
applied to detect topological phase transitions in ID and 
2D spin systems 26,28,29 , as well as in a noninteracting 
fermionic system 26 . 

In our QMC simulations, to detect the correlated QSH 
phase and the interaction-driven phase transition in the 
KMH model, we prepare |fi) as the wave function of Eq. 1 
with U = 0 but keep td different from t. At the nonin¬ 
teracting level, with finite A, td/t will drive a topological 
phase transition between QSH and trivial band insula¬ 
tor at td = 2f 37,39,40 ; therefore throughout this paper we 
choose |fi) with A = 0.2 1 and td = lOOt, which guarantees 
it is a topologically trivial band insulator. On the other 
hand, |T) is prepared as the ground-state wave function 
of the interacting Hamiltonian in Eq. 1 with td = t. In 
the quantum Monte Carlo simulation, it is prepared as 
|T) = e~ eH |Tr), where |Tt) is the wave function of non¬ 
interacting Hamiltonian in Eq. 1 with U = 0, A = 0.2t 
and td = t. The projection operator e~ eH is applied onto 
(T-t) in quantum Monte Carlo sampling such that when 
the projection parameter O is sufficiently large, the QMC 
ensemble average guarantees | SR) is the ground state of 


the interacting Hamiltonian H. In most of the simula¬ 
tions, we set 0 = 50f. 



FIG. 3. (Color online) 1/|C 4 AS | as a function of td in |T). 
The linear divergence of |Ck AS | around the M point holds 
robust until td > 2 1. We can use the divergent to nondivergent 
behavior of |C,( AS | to determine the critical point precisely in 
this noninteracting case. 


In this paper, we define the strange correlator in the 
momentum space. The strange correlator in the single¬ 
particle channel for spin flavor cr is then defined as 


s~iO _ 

°k AB ~ 


(M c kAa- Ck - Bcr l^ < ) 

m 


( 3 ) 


where c' kA(J = j- JA e zk ' Ri ’ A cj A a with k inside the Bril- 
louin zone (BZ) shown in Fig. 2, and A, B are the two 
sublattices of the honeycomb lattice in Fig. 1 (a). The 
schematic plot of Fig. 1 (b) depicts the idea of the strange 
correlator in the KMH model, on the left-hand side, and 
the wave function |fi) is a trivial band insulator (with 
spin Chern number C s = 0); on the right-hand side, the 
projection operator e~ eH guarantees |T) = e _e77 |Tr) is 
the many-body ground state wave function of the KMH 
Hamiltonian at certain U/t , although the trial wave func¬ 
tion |\f r j’) is noninteracting (with spin Chern number 
C s = 1). In this way, as we gradually increase the in¬ 
teraction strength U/t in the KMH Hamiltonian, the na¬ 
ture of |T) will change from QSH at weak interaction 
(U < U c ) to AFMI at strong interaction (U > U c ). 

We also measure the strange correlator in the spin and 
Cooper-pair channels, respectively, as follows: 


c± 

‘- , kA J 4 — 


D^AA — 


WZaS*a\V) 

<fi|T> ’ 

(filA^A^lT) 

m 


( 4 ) 

( 5 ) 


These are two-particle strange correlators in particle- 
hole and particle-particle channels, respectively, where 
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FIG. 4. (Color online) (a). Finite size scaling of the xy anti¬ 
ferromagnetic structure factor for various values of U/t, with 
linear system size L going to 27. The extrapolated values of 
magnetic moment m xy are plotted in (b). 


^ and A kA 


' L e i -*i. 




with S+ A = cj A +Ci } A,i flipping spin in sublattice A of 

unit cell i, and A| a = c] A A ^ creating a Cooper pair 
of a spin singlet in sublattice A of unit cell i. 

Although the magnetic nature of the QSH-to-AFMI 
transition has been studied thoroughly 35,40 , here we find 
the topological nature of this transition is well captured 
by the strange correlators in single- and two-particle sec¬ 
tors. As will be explained later, the QMC computation 
of the strange correlator is more efficient and robust than 
the QMC simulations with either open boundary condi¬ 
tions (OBCs) to directly probe the edge modes 33,34 , or 
measurements of the entanglement spectrum, where one 
has to bifurcate the already small finite-size system and 
analytically continue the imaginary-time data 22,23 . 


III. NUMERICAL RESULTS AND 
DISCUSSIONS 

A. Single-particle strange correlator 

We first apply the single-particle strange correlator to 
detect the topological phase transition driven by td at 
the noninteracting limit. In Fig. 3, we set U = 0 but 
gradually increase td in |*F). One clearly sees that when 
td < 2t, the strange correlator is linearly diver¬ 

gent at one M point in the Brillouin zone, which is consis¬ 
tent with the prediction in Ref. 26. When td > 2 1, both 
= lOOt)| and > 2 1) become a topological triv¬ 
ial state ( C s = 0), the divergence of is removed. 

VICkAsI shows an upturn behavior around kj/. 

Before we move on to the strange correlator in the 
interacting case, we first look at the phase transition from 
QSH to AFMI from the magnetic perspective. Figure. 4 
(a) shows the 1 /L extrapolation of the antiferromagnetic 
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FIG. 5. (Color online) The contour plot of single-particle 
strange correlator |C^ AS | with increasing Hubbard interac¬ 
tion U/t. The finite system size used here is L = 21. The 
k-space area in the four panels is the same as the dashed 
region in the inset of Fig. 2, which is the whole BZ. 


structure factor: 

S af = ^EE (StaSj;* + S~ a S+ a ) (6) 

(i,j) ql=A,B 

for various values of U/t. (• • •) indicates the QMC av¬ 
erage with |'F) on both sides of the observable; hence 
S X AF is not measured as strange correlator but as a reg¬ 
ular QMC correlator. From the extrapolated values of 
L —> oo, one can see the xy antiferromagnetic order sets 
in at around U c « 5.7 1 which is consistent with previous 
QMC results 36,42 . The corresponding magnetic moment 
is obtained as m xy = yjS^p/L 2 , and its value is plotted 
as a function of U/t in Fig. 4 (b). The appearance of 
magnetic long-range order breaks the time-reversal sym¬ 
metry and destroys the bulk topological state. Previous 
theoretical and numerical studies show that the counter 
propagating edge modes in the QSH phase are expected 
to become gapped exactly at the point where long-range 
magnetic order in the bulk breaks time-reversal symme¬ 
try 48 . 

After determining the critical U c from a magnetic per¬ 
spective, we monitor the single-particle strange correla¬ 
tor \C^ An \ as a function of U/t in the whole Brillouin 
zone. For a global view, Fig. 5 shows the contour plot 
of the strange correlator \C^ An \ with increasing U/t for 
a fixed system size L = 21. When U is small, there is a 
clear singularity at one and only one of the time-reversal- 
symmetric M points. In the thermodynamic limit, the 
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single-particle strange correlator is still divergent at one 
M point in Figs. 5 (a-c). When U > U c [see Fig. 5 (d)], 
there will be no divergence in the single-particle strange 
correlator. 

A careful analysis of 1/|C^ AB | along the high- 
symmetry path K—»M—is shown in Fig. 6. In Figs. 6 
(a-d), the single-particle strange correlator still shows 
divergent tendency at the M point with the finite sys¬ 
tem size up to L = 27. When U > U c [see Figs. 6 (e) 
and 6(f)], where the detected target state |'F) becomes 
topologically trivial, we clearly see the upturn behavior 
around the M point in 1/|C^ AB |. 

To give a better understanding of the results in Fig. 6, 
we turn to the helical Luttinger liquid theory. Based on 
the space-time rotation interpretation of the strange cor¬ 
relator, we can analyze the single-particle strange cor¬ 
relator using the helical Luttinger liquid theory at the 
(1+1)D boundary 49-52 , according to which the real-space 
strange correlator in the single-particle sector scales as 

C a tAB ~ r- 9 / 2 - 1 / 2 *, (7) 

where g is the Luttinger parameter related to U/t, g G 
[0,1]. After Fourier transform to the momentum space, 
it becomes 


where k = |k — kM|- Unlike the noninteracting case, the 
single-particle strange correlator in the momentum space 
may actually stop diverging before the QSH to AFMI 
transition point. To see this point more clearly, the criti¬ 
cal g c can be solved from the equation g/2 + l/2g—2 = 0, 
which gives g c = 2 — \/3 ~ 0.268. If g < g c , there would 
be no divergent behavior around the M point in the mo¬ 
mentum space of the single-particle strange correlator, 
although the real-space strange correlator still obeys a 
power-law decay. For g > g c , the power-law divergent 
behavior of the single-particle strange correlator around 
the M point clearly signifies that the interacting QSH 
phase and the trivial band insulator belongs to distinct 
SPT phases, and the two states must be separated by 
gapless fermion edge modes when they are adjacent in 
the space. From the data in Fig. 6, the divergent be¬ 
havior persists up to U = 5.5 1, which is very close to 
the quantum critical point extracted from previous QMC 
simulations. From Fig. 6 here and Fig. 8 in Sec. IIIB, 
we can see that the divergent exponent of the single¬ 
particle strange correlator is reduced by the interaction, 
which cannot be captured by the noninteracting topolog¬ 
ical phase transition in Fig. 3 and is clearly beyond the 
mean-field level. 


C£ab ~ k^ 2+1 ^~ 2 , ( 8 ) 


□ L=9 o L=12 a 1=15 o L=18 o L=21 * L=27 



FIG. 6. (Color online) The inverse amplitude of single-particle strange correlator l/|Cj[ AS | along the high-symmetry path for 
various U/t and system sizes. When the interaction U/t < 5.5 1, see (a-d), there is a divergent tendency in |Ck AS | around the M 
point. However, the divergent exponent is reduced due to the correlation effects according to the helical Luttinger liquid theory 
(see the main text). In the AFMI regime (e-f), l/|Cj[ AB | shows upturn behavior around the M point; there is no divergence at 
all. 


We notice that the data points exactly at k = k m in Figs. 6 (a-d) suddenly jump up and have larger error 
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bars. This is unphysical, and we will discuss the behavior 
of C^ AB in the presence of small antiferromagnetic (AF) 
order Agpw around the M point in a mean-field context 
in Appendix A, where this unphysical singularity at k = 
k m will be understood. 

We want to stress that based on the Luttinger liquid 
theory the single-particle strange correlator, and equiva¬ 
lently, the single particle Green’s function at the physi¬ 
cal edge of the system, always follow a power-law decay 
before the system develops a true long-range order in 
the bulk. This is mainly because when the bulk is fully 
gapped, all the low-energy physics occur at the boundary 
of the system. Then, based on the Mermin-Wagner the¬ 
orem 53 , continuous symmetries cannot be spontaneously 
broken in a (1 + 1)D system, and without a true long- 
range correlation of magnetic or superconductor order 
parameter, the fermions at the boundary should remain 
gapless (though still strongly interacting). 

The technical advantage of the strange correlator in 
QMC over other numerical diagnoses of interacting TIs 
is manifestly presented, i.e., we have performed simula¬ 
tions on finite-size systems with periodic boundary con¬ 
ditions (PBCs) yet still are able to extract information 


on the edge modes, which, in the past, could only be ob¬ 
tained with systems with OBCs 33,34 . It is well known 
that QMC simulations with OBCs suffer from greater 
finite-size effects, but apparently the strange correlator 
avoids this difficulty. Moreover, direct probe of edge 
modes with OBCs requires an analytical continuation of 
the imaginary-time Green’s function, i.e., from G(k, r) to 
A(k,w), and that usually renders ambiguity in the real- 
frequency data. However, with the strange correlator, 
we need to measure only the static (equal time) single¬ 
particle Green’s function in the PBC system, which is the 
easiest and most reliable observables in the QMC simula¬ 
tions. Third, as mentioned in the Introduction, in com¬ 
parison with measurements of the entanglement spec¬ 
trum to detect the interaction-driven topological tran¬ 
sition 21-23 , the strange correlator is also physically more 
transparent and technically more robust, as in the entan¬ 
glement spectrum measurements one has to bifurcate the 
already small finite-size system and analytically continue 
the imaginary-time data, whereas in the strange correla¬ 
tor both problems are avoided. Hence, at the technical 
level, to the best of our knowledge, the strange corre¬ 
lator is indeed the easiest diagnosis of the topological 
states and the topological quantum phase transition in 
interacting systems. 
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FIG. 7. (Color online) (a,c) Spin and pairing strange correlators |Sr aa I and l-DrA^I as a function of system size L for various 
U/t. The inset of (a) is a zoomin at small U/t, with a logarithmic fit (magenta solid line) of the data at U = 0 according to 
Eq. 13. The dark yellow and red solid lines in the main panel of (a) are power-law fits according to Eq. 15 at U = 5.75f(~ U c ) 
and U = 6 t. The logarithmic fit (magenta solid line) in (c) also follows Eq. 14 at U = 0. (b,d) Spin and pairing regular 
correlation functions ('I'|S))"S ; 7|4*} and A).Afc|4/) as a function of L for various U/t. The dark yellow and red solid lines in 
(b) are power-law fits according to Eq. 15 at U = 5.75t(~ U c ) and U = 6f. The inset of (b) is a zoomin of the regular spin 
structure factor at small U/t, showing that they are independent of L when U < U c . 


B. Two-particle strange correlators 

In this section, we discuss the QMC results on strange 
correlators in the two-particle sector, i.e., the spin and 
pairing strange correlators in the presence of interaction. 

Again, based on the space-time rotation interpreta¬ 
tion of the strange correlator, we can likewise analyze 


the spin and pairing strange correlators using the helical 
Luttinger liquid theory at the (1 + 1)D boundary 48-51 . 
According to the theory, the real-space strange correla¬ 
tor in the spin and pairing sectors scale as 

(9) 

D rAA ~ v- 2, \ (10) 






















where g is the Luttinger parameter. After Fourier trans¬ 
form to the momentum space, they become 

^AA~k 2 ®- 2 , (~ L 2 ~ 2a at k = T), (11) 

AcAA~ k 2/ff - 2 , (~ L 2_2/s at k = r). (12) 

where k = |k — kp|. In the noninteracting limit (U = 0) 
<7 = 1, and as we increase U/t towards U c , g will become 
smaller and smaller, eventually vanish at the transition 
point. 

To better understand the behavior in each limit, let us 
start with <7 = 1 (U = 0), and we have 

SkAA ~ k° ~ ln(k), (~ ln(L) at k = T), (13) 
D^ aa ~ k° ~ ln(k), (~ ln(L) at k = T). (14) 

Such a logarithmic growth in L fits our calculated data 
in Figs. 7 (a) and 7 (c) for the U = 0 cases very well. 
The logarithmic growth is in strong contrast to the reg¬ 
ular spin and pairing correlators, as shown in Figs. 7 (b) 
and 7 (d), which, at U = 0, are independent of system 
size L , meaning both spin and pairing correlations are 
exponentially short ranged in real space, corresponding 
to the QSH insulator with a bulk gap. 

On the other hand, near the QSH-to-AFMI transition 
point, g = 0 (U ~ U c ), we have 

'-’kAA ~ ^ i (~ L at k = r), (15) 

A<aa ~ k°°, (~ L~°° ~ e~ L at k = F). (16) 

As we can see, the calculated data in Fig. 7 (a) at U = 6t 
indeed diverge as L 2 in the thermodynamic limit. More 
interestingly, such a divergence is the same as the one 
shown by the regular spin correlator inside the AFMI 
phase, as shown in Fig. 7 (b) at U = 6 1. This is because 
when U = 6 1, the ground state wavefunction |’F) in Eq. 4 
is already in the AFMI phase; the spin strange correlator 
is then similar to the spin regular correlator, because 
both of them pick up the long-range spin-spin correlation. 
In Fig. 7 (c) at U = 6 t, the pairing strange correlator 
decays exponentially to a constant, also consists with the 
prediction in Eq. 16. 

Between the limits of g = 1 (U = 0) and <7 = 0 
(U = f/ c ), we can fit the spin and pairing strange cor¬ 
relator data in Figs. 7 (a) and 7 (c) with the Luttinger 
liquid theory prediction in Eqs. 11 and 12, to extract 
the Luttinger parameter g. The extracted g values as a 
function of U/t are shown in Fig. 8. One can see that 
g continuously decreases from 1 to 0, which accounts for 
the increasing electron-electron correlation. The dashed 
line in Fig. 8 highlights the g c smaller than which the 
single-particle strange correlator stops diverging, as dis¬ 
cussed in Sec. Ill A. 

IV. SUMMARY AND OUTLOOK 

In summary, we have employed large-scale QMC sim¬ 
ulations to study the single-and two-particle strange cor¬ 
relators in a realistic model for interacting topological 



FIG. 8. (color online) The Luttinger parameter g extracted 
from the spin strange correlator in Fig. 7 (a) following Eq. 11. 
Below the critical value g c , the single-particle strange correla¬ 
tor in the momentum space is no longer divergent at M point. 


insulators. We demonstrate that the interaction-driven 
topological-to-trivial quantum phase transition can be 
well captured by the strange correlators. Although larger 
system sizes might be needed for detailed information 
very close to the critical point, our results show that the 
strange correlator is a powerful and promising tool to 
diagnose the topological insulator with interaction. 

The technical advantages of the strange correlator in 
numerical studies (especially QMC simulations) on in¬ 
teracting fermionic and bosonic SPT states are obvious. 
As one needs to measure only static correlations in the 
bulk system, there is no need to apply OBCs to actually 
probe the spatial edges, no need to apply analytical con¬ 
tinuation to access real-frequency data, and no need to 
bifurcate the already small finite-size systems for entan¬ 
glement measurements. In short, the strange correlator 
is much easier to implement and robust in practical nu¬ 
merical performance. 

As for future applications, the QSH insulator discussed 
in our work has a full spin S z conservation, which has a 
Z classification instead of a Z 2 classification for the cases 
with time-reversal symmetry but no S z conservation. In 
Ref. 26 the strange correlator was tested for a noninter¬ 
acting QSH insulator with a sizable Rashba spin-orbit 
coupling, which does have a Z 2 classification. We expect 
that the same strange correlator is still applicable to the 
interacting QSH insulator with Rashba spin-orbit cou¬ 
pling as well, except now that the two electron operators 
in the strange correlator equation. 3 do not have to have 
the same spin, since the spin conservation is broken by 
the Rashba term. 

As we mentioned in the Introduction, in all dimensions 
interaction can change or reduce the classification of some 
topological insulators, for example, interaction may triv¬ 
ialize some topological insulators that are nontrivial in 
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the noninteracting limit. This means that in this case 
the strange correlator should be power-law or long range 
correlated without interaction, but becomes short-ranged 
due to interaction, possibly even without going through 
any bulk phase transition. We will leave this to future 
study. 
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Appendix A: A Mean-Field Calculation of the 
Strange Correlator 

To facilitate the understanding of the behavior of the 
single-particle strange correlator in the QSH insulator 
to xy AFMI transition in the KMH model, below, we 
also provide a mean-field-level calculation of 1 /\C^ AB \ 
by introducing the SDW order parameter, Asdw- The 
mean-field Hamiltonian can be written as 

HmF = ~ ^ ' tij C lcr C jcr + ^ ' v ij c iot a a/3 C j0 

«*?» 

- A S dw ( _ ) z 4a (T S/3 c i/3> 

i,a ,/3 

where Asdw is the SDW gap. Here we set A = 0.22. 
If t,i = t and Asdw = 0, L/mf describes the QSH in¬ 
sulator. The trivial band insulator can be obtained by 
tuning td > 2 1. The strong interacting AFMI can be 
phenomenologically modeled by a finite Asdw term in 
the mean-field theory, which breaks the spin C7(l) sym¬ 
metry and describes the spin ordered antiferromagnetic 
state. Here we assume the xy spin order lies in the spin-a; 
direction. 

We calculated the inverse strange correlator ^/\Cl. AB \ 
with the state IT) tuned by the mean-field parameter 
Asdw- The result is shown in Fig. 9. When Asdw = 0, 
|T) is the QSH state, the inverse strange correlator 
Vl^kAsI ~ |k — k M | follows the linear behavior around 
the M point, which implies the power law behavior of 
the strange correlator ~ |k — kA/| - “ with a = 1. 

However, beyond the mean-field theory, the interaction 
can modify this power a, so the strange correlator can de¬ 
viates from the a = 1 behavior in the momentum space, 



FIG. 9. The inverse strange correlator 1/|Gj( ab | along the 
high-symmetry path in the mean-held theory. The state |T) 
is replaced by an SDW insulator controlled by Asdw- 


as shown in Eq. 8 and Fig. 8 in the main text. But the 
power-law behavior of the strange correlator in the real 
space is still expected to survive in the whole QSH phase. 

As Asdw is turned on, l/|Ck AB | will be lifted from 
zero at the M point and replaced by a small peak. The 
stronger SDW order will leads to the earlier upturn of 
the curve as approaching to the M point. The upturn 
behavior around the M point can be described by 


t _ (k + »mn)(Ag DW mf! + k(k - imy)(mn + my)) 
kAB 2(Ag DW mQ + k 2 (mn + my) 2 ) 

(A2) 

where k = r>p-|k — k/\/| is the small momentum deviation 
from the M point, and ran and my are, respectively, 
the single-particle mass gaps in the trivial state |D) and 
the QSH state |T). Equation. A2 is derived by a small 
momentum expansion around the M point. As can be 
seen from the denominator, the power law divergence of 
the strange correlator |C^ AB | (as k -A 0) will be removed 
once the SDW order A S dw sets in at the topological 
transition to the AFMI phase. Additionally, according to 
Eq. A2, the limit of Asdw -a 0 and the limit of k — > k m 
do not commute: 


lim Ck AB = 

fc-> o 


i 

2 ’ 


^ _ {k + imn){k - imn) 

lim C k a r — 7 / \ 

a S dw->-o k(mn + my) 


1 

k’ 


(A3) 


If one takes Asdw -a 0 first, then the strange correla¬ 
tor C^ AB indeed follows the 1/k power-law behavior as 
expected on the mean-field level. However, in our QMC 
simulation, we take k —> 0 first due to the presence of 
the AF fluctuations as a result of the finite-size effect, 
so the strange correlator C^ AB approaches another limit 
*/2, which is not divergent. It is this non-commutative 











10 


limit that makes | C^ A B | ill-defined at the M point and 
the data right at the M point meaningless in Fig. 3 (a- 
d). Only when the interaction becomes sufficiently strong 


(the Luttinger parameter g < g c ), the single-particle 
strange correlator will no longer diverge, does the data 
of l/|CkAsl at k = M becomes meaningful. 
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